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Introduction 



LanHEP has been designed as a part of the CompHEP package Jlj], worked out for automatic calculations in 
high energy physics. CompHEP allows symbolic computation of the matrix element squared of any process with 
up to 6 incoming and outgoing particles for a given physical model (i.e. a model defined by a set of Feynman 
rules as a table of vertices in the momentum representation) and then numerical calculation of cross-sections 
and various distributions. 

Main purpose of the new option given by LanHEP program is designing of a new physical model. LanHEP 
program makes possible the generation of Feynman rules for propagators and vertices in momentum representa- 
tion starting from the Lagrangian defined by a user in some simple format very similar to canonical coordinate 
representation. User should prepare a text file with description of all Lagrangian terms in the format close to 
the form used in standard publications. Of course, user has to describe all particles and parameters appearing 
in Lagrangian terms. 

The main LanHEP features are: 

• LanHEP expands expression and combines similar terms; 

• it performs the Fourier transformation by replacing derivatives with momenta of particles; 

• it writes Feynman rules in the form of four tables in CompHEP format as well as tables in LaTeX format; 

• user can define the substitution rules, for example for covariant derivative; 

• it is possible to define multiplets, and (if necessary) their components; 

• user can write Lagrangian terms with Lorentz and multiplet indices explicitly or omit indices (all or some 
of them); 

• LanHEP performs explicit summation over the indices in Lagrangian terms, if the corresponding compo- 
nents for multiplets and matrices are introduced; 

• it allows the user to introduce vertices with 4 fermions and 4 colored particles (such vertices can't be 
introduced directly in CompHEP) by means of auxiliary field with constant propagator; 

• it also can check whether the set of introduced vertices satisfies the electric charge conservation law. 

1 QED by means of LanHEP 

We start from a simple exercises, illustrating the main ideas and features of LanHEP. The first physical model 
is Quantum Electrodynamics. 



model QED/1 . 

parameter ee=0 . 31333 :' elementary electric charge', 
spinor el/El : (electron, mass me=0 . 000511) . 
vector A/A: (photon) . 

let F~mu~mi=deriv~mi*A~mu-deriv~mu*A~nu . 

lterm -1/4* (F'mu'nu) **2 - 1/2* (deriv~mu*A~mu) **2 . 

lterm El* (i*gamma*deriv+me) *el . 

lterm ee*El*gamma*A*el . 



Figure 1: LanHEP input file for the generation of QED Feynman rules 
QED Lagrangian is 

Cqed = --F^F^ + e-f"^ + g e A^e - mee 
and the gauge fixing term in Feynman gauge has the form 

C GF = - l -{d,A»)\ 

Here e(x) is the spinor electron-positron field, m is the electron mass, A^x) is the vector photon field, F^ v = 
d v A^ — d^A u , and g e is the elementary electric charge. 
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The LanHEP input file to generate the Feynman rules for QED is shown in Fig. 

First of all, the input file consists of statements. Each statement begins with one of the reserved keywords 
and ends by the full-stop '.' symbol. 

First line says that this is a model with the name QED and number 1. This information is supplied for 
CompHEP, the name QED will be displayed in its list of models. In CompHEP package each model is described 
by four files: 'varsN.mdl', 'funcN.mdl', 'prtclsN.mdl', TgrngnN.mdP , where N is the very number specified in 
the model statement. 

The model statement stands first in the input file. If this statement is absent, LanHEP does not generate 
four standard CompHEP files, just builds the model and prints diagnostic if errors are found. 

Second line in the input file contains declaration of the model parameter, denoting elementary electric charge 
g e as ee. For each parameter used in the model one should declare its numeric value and optional comment (it 
is also used in CompHEP menus). 

The next two lines declare particles. Statement names spinor , vector correspond to the particle spin. So, 
we declare electron denoted by el (the corresponding antiparticle name is El) and photon denoted by A (with 
antiparticle name being A, since the antiparticle for photon is identical to particle). 

After particle name we give in brackets some options. The first one is full name of the particle, used in 
CompHEP; the second option declares the mass of this particle. 

The let statement in the next line declares the substitution rule for symbol F, which will be replaced in the 
further Lagrangian terms by the expression given in this statement. 

Predefinite name deriv, reserved for the derivation will be replaced after the Fourier transformation by 
the momentum of the particle multiplied by —i. 

The rest of the lines describe Lagrangian terms. Here the reserved name gamma denotes Dirac's 7-matrices. 

One can see that the indices are written separated with the caret symbol Note that in the last two lines 
we have omitted indices. It means that LanHEP restores omitted indices automatically. Really, one can type 
the last term in the full format: 

lterm ee*El~a*gamma~a~b~mu*A~mu*el~b. 

It corresponds to g e e a 7^ b e(,j4 M with all indices written. Note that the order of objects in the monomial is 
important to restore indices automatically. 



2 QCD 

Now let us consider the case of the Quantum Chromodynamics. The Lagrangian for gluon fields reads 

Lym = -\F a ^ v F^ 

where 

r III* ~ °^ (j v ~~ ~ 9s J 

G®(x) is the gluon field, g s is a strong charge and f abc are purely imaginary structure constants of SU(3) color 
group. 

The quark kinetic term and its interaction with the gluon has the form 

L f = q% 7 M <9 M q t + g s A? q % 7 M q 3 ■ G£ , 

where A" are Gell-Mann matrices. 

Gauge fixing terms in Feynman gauge together with the corresponding Faddev-Popov ghost term are 

where (c, c) are unphysical ghost fields. 

The corresponding LanHEP input file is shown in Fig. |^. 

Since QCD uses objects with color indices, one has to declare the indices of these objects. There are three 
types of color indices supported by LanHEP. These types are referred as color c3 (color triplets), color c3b 
(color antitriplets), and color c8 (color octets). One can see that color c3 index type appears among the 
options in the quark q declaration, and the color c8 one in the gluon G declaration. Antiquark Q has got color 
index of type color c3b as antiparticle to quark. LanHEP allows convolution of an index of type color c3 
only with another index of type color c3b, and two indices of type color c8. Of course, in Lagrangian terms 
each index has to be convoluted with its partner, since Lagrangian has to be scalar. 
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model QCD/2. 

parameter gg=l . 117 :' Strong coupling'. 

spinor q/Q: (quark, mass mq=0.01, color c3) . 

vector G/G:(gluon, color c8, gauge). 

let F~mu~nu~a = deriv~nu*G~mu~a - deriv~mu*G~nu~a - 

gg*f _SU3~a~b~c*G~mu~b*G~nu~c . 
lterm -F**2/4- (deriv*G) **2/2 . 
lterm Q*(i*gamma*deriv+mq)*q. 
lterm i*gg*f _SU3*ccghost (G) *G*deriv*ghost (G) . 
lterm gg*Q*gamma*lambda*G*q. 



Figure 2: Input file for the generation of QCD Feynman rules 



Table 1: QCD Feynman rules generated by LanHEP in LaTeX output format 



Fields in the vertex 


Variational derivative of Lagrangian by fields 


G>p G.Cq G.c r 

Qap Qbq G pr 
G pp Gyq G pr 

G p,p Gyq G p r G(j S 


-gg-Psfpqr 

gg ■ iab x pq 

ggfpqr (p^ p - P p 2 g pv - v%9 vp + p p g pv + p^g vp - pW ip ) 
gg 2 (g pp g ua fpqtfrst - g^g" p f Pq tfrst + g» v g pa fprtUt 

+g PV 9 Pa fpstf q rt - g^g" P fprtfqst - g^g^fpstfqrt) 



LanHEP allows also to use in the Lagrangian terms a predefined symbol lambda with the three indices of 
types color c3, color c3b, color c8 corresponding to Gell-Mann A-matrices. Symbol f_SU3 denotes the 
structure constant f abc of color SU(3) group (all three indices have the type color c8). 

Option gauge in the declaration of G allows to use names ghost (G) and ccghost(G) for the ghost fields c 
and c in Lagrangian terms and in let statements. 

Table 1 shows Feynman rules generated by LanHEP in LaTeX format after processing the input file presented 
in Fig. H Four gluon vertex rule is indicated in the last line. Note that the output in CompHEP format has no 
4-gluon vertex explicitly; it is expressed effectively through 3-leg vertices by the constant propagator of some 
auxiliary field (see Appendix [A| for more details). 

3 Syntax 

The LanHEP input file is the sequence of statements, each starts with a special identifier (such as parameter , 
lterm etc) and ends with the full-stop '.' symbol. Statement can occupy several lines in the input file. 
This section is aimed to clarify the syntax of LanHEP input files, i.e. the structure of the statements. 

3.1 Constants and identifiers 

First of all, each word in any statement is either an identifier or a constant. 

Indentfiers are the names of particles, parameters etc. Examples of identifiers from the previous section are 
particle names 

el El A q Q G 

The first word in each statement is also an identifier, defining the function which this statement performs. The 
identifiers are usually combinations of letters and digits starting with a letter. If an identifier doesn't respect 
this rule, it should be quoted. For example, the names of bosons must be written as 'W+' and 'W-', since 
they contain '+' and '-' symbols. 
Constants can be classified as 

• integers: they consist of optional sign followed by one or more decimal digits, such as 
1 -1 123 -98765 
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Integers can appear in Lagrangian terms, parameter definition and in other expressions. 

• Floating point numbers include optional sign, several decimal digits of mantissa with an embedded period 
(decimal point) with at least one digit before and after the period, and optional exponent. The exponent, 
if present, consist of letter E or e followed by an optional sign and one or more decimal digits. The valid 
examples of floating point numbers are 

1.0 -1.0 0.000511 5.11e-4 

Floating point numbers are used only as parameter values (coupling constants, particle masses etc). They 
can not be explicitly used in Lagrangian terms. 

• String constants may include arbitrary symbols. They are used as comments in parameter statements, 
full particle names in the declaration of a particle, etc. Examples from the previous section are 

electron 
photon 

If a string constant contains any character besides letters and digits or doesn't begin with a letter, it 
should be quoted. For example, the comments in QED and QCD input files (see previous section) contain 
blank spaces, so they are quoted: 

'elementary electric charge' 
'Strong coupling' 

3.2 Comments 

User can include comments into the LanHEP input file in two ways. First, symbol '°/ ' denotes the comment 
till the end of current line. Second way allows one to comment any number of lines by putting a part of input 
file between '/*' (begin of comment) and '*/' (end of comment) symbols. 

4 Objects in the expressions for Lagrangian terms 

Each symbol which may appear in algebraic expressions (names of parameters, fields, etc) has a fixed order of 
indices and their types. If this object is used in any expression, one should write its indices in the same order 
as they were defined when the object has been declared. 

Besides the indices types corresponding to color SU(3) group: color c3 (color triplet), color c3b (color 
antitriplet) and color c8 (color octet) described in the previous example, there are default types of indices 
for Lorentz group: vector, spinor and cspinor (antispinor). User can also declare new types of indices 
corresponding to the symmetries other than color SU(3) group. In this case any object (say, particle) may have 
indices related to this new group. This possibility will be described in Section 0. 

If an index appears twice in some monomial of an expression, LanHEP assumes summation over this index. 
Types of such indices must allow the convolution, i.e. they should be one of the pairs: spinor and cspinor, 
two vector, color c3 and color c3b, two color c8. 

In general the following objects are available to appear in the expressions for a Lagrangian: integers and 
identifiers of parameters, particles, specials, let-substitutions and arrays. 

There are also predefined symbols i, denoting imaginary unit i (i 2 — — 1) and Sqrt2, which is a parameter 
with value the \/2. 

4.1 Parameters 

Parameters are scalar objects (i.e. they have no indices). Parameters denote coupling constants, masses and 
widths of particles, etc. To introduce a new parameter one should use the parameter statement, which has the 
generic form 

parameter name=value: comment. 

• name is an identifier of newly created parameter. 

• value is an integer or floating point number or an expression. One can use previously declared parameters 
and integers joined by standard arithmetical operators ' + ', '-', '*', '/', and ' ** ' (power). 
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• comment is an optional comment to clarify the meaning of parameter, it is used in CompHEP help 
windows. Comment has to be a string constant, so if it contains blank spaces or other special characters, 
it must be quoted (see Section 3). 

4.2 Particles 

Particles are objects to denote physical particles. They may have indices. It is possible to use three statements 
to declare a new particle, at the same time the second and the third statements define the corresponding Lorentz 
index: 

scalar P/aP: (options) . 
spinor P/aP: (options) . 
vector P/aP: (options) . 

P and aP are identifiers of particle and antiparticle. In the case of truly neutral particles (when antiparticle is 
identical to the particle itself) one should use the form P/P with identical names for particle and antiparticle. 
It is possible to write only the particle name, e.g. 

scalar P: (options) . 

In this case the name of corresponding antiparticle is generated automatically. It satisfies the usual CompHEP 
convention, when the name of antiparticle differs from particle by altering the case of the first letter. So for 
electron name el automatically generated antiparticle name will be El. If the name contains symbol '+' it is 
replaced by '-' and vice versa. 

The option is comma-separated list of options for a declared particle, and it may include the following items: 

• the first element in this list must be the full name of the particle, (e.g. electron and photon in our 
example.) Full name is string constant, so it should be quoted if it contains blanks, etc. 

• mass param=value defines the mass of the particle. Here param is an identifier of new parameter, which is 
used to denote the mass; value is its value, it has the same syntax as in the parameter statement, comment 
for this new parameter being generated automatically. If this option omitted, the mass is assumed to be 
zero. 

• width param^value declares the width of the particle. It has the same syntax as for mass option. 

• itype is a type of index of some symmetry; one can use default index types for color SU (3) group (see 
QCD example in Section 2). It is possible to use user-defined index types (see Section 0) and Lorentz 
group indices vector, spinor, cspinor. 

• left or right say that the massless spinor particle is an eigenstate of (1 — 7 5 )/2 or (1 +7 5 )/2 projectors, 
so this fermion is left-handed or right-handed one. 

• gauge declares the vector particle as a gauge boson. This option generates corresponding ghosts and 
goldstone bosons names for the named particle (see below). 

When a particle name is used in any expression (in Lagrangian terms), one should remember that the first 
index is either vector or spinor one (of course, if this particle is not a Lorentz scalar). Then the indices follow 
in the same order as index types in the options list. So, in the case of quark declaration (see QCD example) 
the first index is spinor, and the second one is color triplet. 

There are several functions taking particle name as an argument which can be used in algebraic expressions. 
These functions are replaced with auxiliary particle names, which are generated automatically. 

• Ghost field names in gauge theories are generated by the functions ghost (name) — > 'namc.c' and 
ccghost (name) — * 'namc.C (see for instance Table 1). Here and below name is the name of the 
corresponding gauge boson. 

• Goldstone boson field name in the t'Hooft-Feynman gauge is generated by the function gsb(name) — > 
' name . f ' . 

• The function anti(name) generates antiparticle name for the particle name. 

• The name for a charge conjugated spinor particle ip c — Ctp T is generated by the function cc(name) — > 
' name . c ' . Charge conjugated fermion has the same indices types and ordering, however index of spinor 
type is replaced by the index type cspinor and vice versa. 
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• vev(expr) is used in Lagrangian terms for vacuum expectation values. Function vev ensures that 
deriv*vev(expr) is zero. In other words, vev function forces LanHEP to treat expr as a scalar par- 
ticle which will be replaced by expr in Feynman rules. 



4.3 Specials 

Besides parameters and particles other indexed such as 7-matrices, group structure constants, etc may appear 
in the Lagrangian terms. We refer such objects as specials. 
Predefined specials of Lorentz group are: 

• gamma stands for the 7-matrices. It has three indices of spinor, cspinor and vector types. 

• gamma5 denotes 7 s matrix. It has two indices of spinor and cspinor types. 

• moment has one index of vector type. At the stage of Feynman rules generation this symbol is replaced 
by the particle moment. 

• deriv is replaced by —ip^, where p^ is the particle moment. It has one vector index. 
Specials of color SU (3) group are: 

• lambda denotes Gell-Mann A-matrices. It has three indices: color c3 , color c3b and color c8. 

• f _SU3 is the SU(3) structure constant. It has three indices of color c8 type. 
Note that for specials the order of indices types is fixed. 

User can declare new specials with the help of a facility to introduce user-defined indices types (see Section 0) . 

4.4 Let-substitutions 

LanHEP allows the user to introduce new symbols and then substitute them in Lagrangian terms by some 
expressions. Substitution has the generic form 

let name=expr. 

where name is the identifier of newly defined object. The expression has the same structure as those in 
Lagrangian terms, however here expression may have free (non-convoluted) indices. 

Typical example of using a substitution rule is a definition of the QED covariant derivative as 

let Deriv~mu=deriv~mu + i*ee*A~mu. 

corresponding to = + igeA^. 

More complicated example is the declaration er^ = 4(7^7^ — ^j^)/2 matrices: 

let sigma~a~b~mu~nu = i* (gamma" a~c~mu*gamma~c~b~nu - gamma" a~c~nu*gamma~c~b~mu) II . 

Note that the order of indices types of new symbol is fixed by the declaration. So, first two indices of sigma 
after this declaration are spinor and antispinor, third and fourth are vector indices. 

4.5 Arrays 

LanHEP allows to define components of indexed objects. In this case, convolution of indices will be performed 
as an explicit sum of products of the corresponding components. 
Object with explicit components has to be written as 

{exprl, expi2 exprN }~i 

where expressions correspond to components. All indices of components (if they present) have to be written at 
each component, and the index numbering components has to be written after closing curly bracket. Of course, 
all the components must have the same types of free (non-convolved) indices. 

Arrays are usually applied for the definition of multiplets and matrices corresponding to broken symmetries. 

Typical example of arrays usage is a declaration of electron-neutrino isospin doublet 11 (and antidoublet 
LI) ' 

let ll~a~I = { nl"a, el~a}~I, Ll~a~I = { Nl~a, El"a}"I. 

Here we suppose that nl was declared as the spinor particle (neutrino), with the antiparticle name Nl. 

Matrices can be represented as arrays which have other arrays as components. However, it is more convenient 



to declare them with omitted indices, see Section 3.4 (the same is correct for multiplets also). 

It is possible also to use arrays directly in the Lagrangian terms, rather than only in the declaration of 
let-substitution. 
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5 Structure of expressions for the Lagrangian 



When all parameters and particles necessary for introduction of physical model are declared, one can enter 
Lagrangian terms with the help of the lterm statement: 

lterm expr. 

Elementary objects of expression are integers, identifiers of parameters, particles, specials, let-substitutions, 
and arrays. 

These elementary objects can be combined by usual arithmetical operators as 

• exprl+expr2 (addition), 

• exprl-expr2 (subtraction), 

• exprl*expr2 (product), 

• exprl/expr2 (fraction; here expr2 must be a product of integers and parameters), 

• exprl**N (Nth power of exprl; N must be integer). 

One can use brackets '(' and ')' to force the precedence of operators. Note, that indices can follow only 
elementary objects symbols, i.e. if Al and A2 were declared as two vector particles then valid expression for 
their sum is Al~mu+A2~mu, rather than (A1+A2) ~mu. 



5.1 Where-substitutions 

More general form of expressions involves where-substitutions: 
expr where subst. 

In the simple form subst is name=repl or several constructions of such kind separated by comma In the 
form of such kind each instance of identifier name in expr is replaced by repl. 

Note that in contrast to let-substitutions, where-substitution doesn't create a new object. LanHEP simply 
replaces name by repl, and then processes the resulting expression. It means in particular that name can not 
have indices, although it can denote an object with indices: 

lterm F**2 where F=deriv~mu*A~nu-deriv~nu*A~mu. 

is equivalent to 

lterm (deriv~mu*A~nu-deriv~mi*A~mu) **2 . 

The substitution rule introduced by the keyword where is active only within the current lterm statement. 

More general form of where-substitution allows to use several name=repl substitution rules separated by 
semicolon In this case expr will be replaced by the sum of expressions; each term in this sum is produced by 
applying one of the substitution rules from semicolon-separated list to the expression expr. This form is useful 
for writing the Lagrangian where many particles have similar interaction. 

For example, ifu,d,s,c,b,t are declared as quark names, their interaction with the gluon may read as 

lterm gg*anti (psi) *gamma*lambda*G*psi where 

psi=u; psi=d; psi=s; psi=c; psi=b; psi=t. 

The equivalent form is 

lterm gg*U*gamma*lambda*G*u + gg*D*gamma*lambda*G*d + 
gg*U*gamma*lambda*G*u + gg*D*gamma*lambda*G*d + 
gg*U*gamma*lambda*G*u + gg*D*gamma*lambda*G*d. 

Where-substitution can also be used in let statement. In this case one should use brackets: 

let lsub=(expr where wsub=exprl). 

Note that in the previous example we have omitted indices. In the next section we shall describe this 
LanHEP option in details. 
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6 Omitting indices 



Physicists usually do not write all possible indices in the Lagrangian terms. LanHEP also allows a user to omit 
indices. This feature can simplify introduction of expressions and make them more readable. Compare two 
possible forms: 

lterm El~a*gamma~a~b~mu*A~mu*el~b . 

corresponding to g e e a (x)^ b e b (x)A^(x), and 

lterm El*gamma~mu*el*A~nm. 

corresponding to g e e(x)-y^e(x)A fl (x)). Furthermore, while physicists usually write vector indices explicitly in 
the formulas, in LanHEP vector indices also can be omitted: 

lterm -i*ee*El*gamma*el*A . 

Generally speaking, when user omits indices in the expressions, LanHEP faces two problems: which indices 
were omitted and how to convolute restored indices. 

6.1 Restoring the omitted indices 

When the indexed object is declared the corresponding set of indices is assumed. Thus, if the quark q is declared 

as 

spinor q:('some quark', color c3) . 

its first index is spinor and the second one belongs to the color c3 type. If both indices are omitted in some 
expression, LanHEP generates them in the correspondence to order (spinor, color c3). However, if only one 
index is written, for example in the form q~a, LanHEP has to recognize whether the index a is of color c3 or 
of spinor types. 

To solve this problem LanHEP looks up the list of indices omitting order. By default this list is set to 

[spinor, color c3, color c8, vector] 

The algorithm to restore omitted indices is the following. First, LanHEP assumes that user has omitted 
indices which belong to the first type (and corresponding antitype) from this list. Continuing the consideration 
of our example with particle q one can see that since this particle is declared having one spinor index (the first 
type in the list) LanHEP checks whether the number of indices declared for this object without spinor index 
equals to the number of indices written explicitly by user. In our example (when user has written q~a) this is 
true. In the following LanHEP concludes that the user omitted spinor index and that explicitly written index 
is of color c3 type. 

In other cases, when the supposition fails if the user has omitted indices of the first type in the list of indices 
omitting order, LanHEP goes to the second step. It assumes that user has omitted indices of first two types 
from this list. If this assumption also fails, LanHEP assumes that user has omitted indices of first three types in 
the list and so on. At each step LanHEP subtracts the number of indices of these types assumed to be omitted 
from the full number of indices declared for the object, and checks whether this number of resting indices equals 
to the number of explicitly written indices. If LanHEP fails when the list of indices omitting order is completed, 
error message is returned by the program. 

Note that if the user would like to omit indices of some type, he must omit all indices of this type (and 
antitype) as well as the indices of all types which precede in the list of indices omitting order. 

For example, if object Y is declared with one spinor, two vector and three color c8 indices, than 

• the form Y~a~b~c~d~e~f means that the user wrote all the indices explicitly; 

• the form Y~a~b~c~d~e means that the user omitted spinor index and wrote vector and color c8 ones; 

• the form Y~a~b means that the user omitted spinor and color c8 indices and wrote only two vector 

ones; 

• the form Y means that the user omitted all indices; 

• all other forms, involving different number of written indices, are incorrect. 
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One can say that indices should be omitted in the direct correspondence with abrupting the list of indices 
omitting order from left to right. 

One could change the list of indices omitting order types by the statement SetDef Index. For example, for 
default setting it looks like 

SetDef Index(spinor, color c3, color c8, vector). 

Each argument in the list is a type of index. 



6.2 Convolution of restored indices 

Omitted index can be convoluted only with some another omitted index. LanHEP expands the expression 
and restore indices in each monomial. LanHEP reads objects in the monomial from the left to the right and 
checks whether restored indices are present. If such index appears LanHEP seeks for the restored index of the 
appropriate type at the next objects. Note, that the program does not check whether the object with the first 
restored index has another restored index of the appropriate type. Thus, if F is declared as let-substitution for 
the strength tensor of electromagnetic field (with two vector indices) then expression F*F (as well as F**2) after 
processing omitted indices turns to implied form F~mu~nu*F~mu~nu rather than F~mu~mu*F~mi~mi. 

This algorithm makes the convolutions to be sensitive to the order of objects in the monomial. Let us look 
again to the QED example. Expression El*gamma*A*el (as well as A*el*gamma*El) leads to correct result where 
vector index of photon is convoluted with the same index of 7-matrix, spinor index of electron is convoluted 
with antispinor index of 7-matrix and antispinor index of positron is convoluted with spinor index of 7-matrix. 
However the expression el*El*A*gamma leads to wrong form el~a*El~a*A~mu*gamma~b~c~mu, because the first 
antispinor index after electron belongs to positron. Spinor indices of gamma stay free (non-convoluted) since no 
more objects with appropriate indices (so, LanHEP will report an error since Lagrangian term is not a scalar). 

Note, that in the vertex with two 7-matrices the situation is more ambiguous. Let's look at the term 
corresponding to the electron anomalous magnetic moment e(x)(7 M 7 iy — j' / -f^)e(x)F fll ,. The correct LanHEP 
expression is 

el* (gamma~mu*gamma~nu - gamma~nu*gamma~mu) *El*F~mu~nu 

Here vector indices can't be omitted, since it lead to the convolution of vector indices of 7-matrices. One can 
see also that the form 

el*El* (gamma~mu*gamma~nu - gamma~nu*gamma~mu) *F~mu~nu 

will correspond to the expression e(x)e(x)Sp('y ti j")F lll/ . Here LanHEP has got scalar Lorentz-invariant expres- 
sion in the Lagrangian term, so it has no reason to report error. 

These examples mean that user should clearly realize how the indices will be restored and convoluted, or he 
has to write all indices explicitly. 



6.3 Let-substitutions 

Other problem arises when the omitted indices stay free, it is a case for the let statement. LanHEP allows only 
two ways to avoid ambiguity in the order of indices types: either user specifies all the indices at the name of new 
symbol and free indices in the corresponding expression, or he should omit all free indices. In the latter case 
the order of indices types is defined by the order of free omitted indices in the first monomial of the expression. 
For example if Al and A2 are vectors and cl and c2 are spinors, the statement 

let d=Al*cl+c2*A2. 

declares new object d with two indices, the first is vector index and the second is spinor index according to their 
order in the monomial Al*cl. Of course, each monomial in the expression must have the same typeset of free 
indices. 



6.4 Arrays 

The usage of arrays with omitted indices allows us to define matrices conveniently. For example, the declaration 
of r-matrices 

' 1 

T ' = l 

can be written as 



7-1 = 1! 1 ' T2 



i \ 








-i ) 


, „-( 


> 


-J 
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let taul = {{0, 1}, {1, 0}}. 
let tau2 = {{0,-i}, {i, 0}}. 
let tau3 = {{1, 0}, {0,-1}}. 

One can see that in such way of declaration a matrix is written " column by column" . 

The declaration of objects with three 'explicit' indices can be done using the objects already defined. For 
example, when r- matrices are defined as before, it is easy to define the vector f = (tx, T2,Ts) as 

let tau = {taul, tau2, tau3} . 

The object tau has three indices, first pair selects the element of the matrix, while the matrix itself is selected 
by the third index, i.e. tau~i~j~a corresponds to r^. 

On the other hand, the declaration of structure constants of a group is more complicated. Declaring such 
an object one should bear in mind that omitting indices implies that in a sequence of components the second 
index of an object is changed after the full cycle of the first index, the third index is changed after the full cycle 
of the second one, etc. For example, a declaration of the antisymmetrical tensor e abc can read as 

let eps = {{{0,0,0}, {0,0,-1}, {0,1,0}}, 
{{0,0,1}, {0,0,0}, {-1,0,0}}, 
{{0,-1,0}, {1,0,0}, {0,0,0}}}. 

One can easily see that the components are listed here in the following order: 

£ 111 £ 211 £ 311 £ 121 £ 221 £ 321 £ 131 £ 233 £ 333 

The declaration of more complex objects such as SU (3) structure constants can be made in the same way. 

7 Declaration of new index types and indexed objects 

7.1 Declaring new groups 

Index type is defined by two keywords^]: group name and representation name. Thus, color triplet index type 
color c3 has group name color and representation name c3. 

LanHEP allows user to introduce new group names by the group statement: 

group gname. 

Here gname is a string constant, which becomes the name of newly declared group. 
Representation names for each group name must be declared by the statement 

repres gname: (rlist) 

where rlist is a comma-separated list of representation names declaration for the already declared group name 
gname. Each such declaration has the form cither mame or rname/ crname. In the first case index which belongs 
to the gname rname type can be convoluted with another index of the same type; in the second case index of 
gname rname type can be convoluted only with an index of gname crname type. 

For example, definition for color SU (3) group with fundamental, conjugated fundamental and adjoint rep- 
resentations looks as: 

group color:SU(3). 
repres color : (c3/c3b,c8) . 

So, three index types can be used: color c3, color c3b, color c8. The convolution of these indices is 
allowed by pairs (color c3, color c3b) and (color c8, color c8) indices. 

7.2 Declaring new specials 

Specials with indices of user-defined types can be declared by means of special statement: 
special name : (islist) . 

Here name is the name of new symbol, and Hist is a comma-separated list of indices types. For example, 
Gell-Mann matrices can be defined as (although color group and its indices types are already defined): 

special lambda: (color c3, color c3b, color c8) . 

To define Dirac's 7-matrices we can use the command 
special gamma: (spinor, cspinor, vector). 

lr The exception is Lorentz group, corresponding indices types are defined by single keyword. 
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7.3 Arrays 

Array, i.e. the object with explicit components, can also have the user-defined type of index. In this case generic 
form of such object is 

{ exprl, expr2 ... ,exprN ; itype } 

where N expressions exprl ... expr N of N components are separated by comma, and itype is an optional index 
type. If itype is omitted LanHEP uses default group name wild and index type wild N, where N is a number 
of components in the array. 

8 Auxiliary statements 

8.1 Orthogonal matrices 

If some parameters appear to be the elements of the orthogonal matrix such as quark mixing Cabbibo-Kobayashi- 
Maskava matrix, one should declare them by the statement 

OrthMatrixC {{an, o 12 , 013} , {021, 022, 023}, {a 3 i, 032, 033}} )■ 

where denote the parameters. Such declaration permits LanHEP to reduce expressions which contain these 
parameters by taking into account the properties of orthogonal matrices. 

Note that this statement has no relation to the arrays; it just declares that these parameters satisfy the 
correspondent relations. Of course, one can declare further a matrix with these parameters as components by 
means of let statement. 

8.2 Including files 

LanHEP allows the user to divide the input file into several files. To include the file Gle, the user should use 
the statement 

read file. 

The standard extension '.mdl' of the file name may be omitted in this statement. 
Another way to include a file is provided by the use statement as 

use Gle. 

The use statement reads the file only once, next appearances of this statement with the same argument do 
nothing. This function prevents multiple reading of the same file. This form can be used mainly to include 
some standard modules, such as declaration of Standard Model particles to be used for writing some extensions 
of this model. 

8.3 Checking electric charge conservation 

LanHEP can check whether the introduced vertices satisfy electric charge conservation law. This option is 
available, if the user declares some parameter to denote elementary electric charge (say, ee in QED example), 
and than indicate, which particle is a photon by the statement 

SetEM (photon, param). 

So, in example of Section 1 this statement could be 

SetEM(A.ee) . 

Electric charge of each particle is determined by analyzing its interaction with the photon. LanHEP checks 
whether the sum of electric charges of particles in each vertex equals zero. 
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9 LaTeX output 



LanHEP generates LaTeX output instead of CompHEP model files if user set -tex in the command line to start 
LanHEP. Three files are produced: 'varsJV.tex', 'prtclsJV.tex' and 'lgrngJV.tex'. The first file contains 
names of parameter used in physical model and their values. The second file describes introduced particles, 
together with propagators derived from introduced vertices. 

The last file lists introduced vertices. LanHEP uses Greece letters fi, v, p... for vector indices, letters a, b, c... 
for spinor ones and p, q,r... for color indices (and for indices of other groups, if they were defined). 

It is possible to inscribe names for particles and parameters to use them in LaTeX output. It can be done 
by the statement 

SetTexName ( [ ident=texname, ... ]). 

Here ident is an identifier of particle or parameter, and texname is string constant containing LaTeX command. 
Note, that for introducing backslash \ in quoted string constant one should type it twice: '\\'. 

For example, if one has declared neutrino with name nl (and name for antineutrino Nl) than the statement 

SetTexName([nl='\\mTe', Nl= ' \\bar{\\nu}~e ' ] ) . 

makes LanHEP to use symbols v e and D e for neutrino and antineutrino in LaTeX tables. 

10 Running LanHEP 

As it was mentioned above, LanHEP can read the model description from the input file prepared by user. To 
start LanHEP write the command 

lhep filename options 

where the possible options are described in the next section. If the filename is omitted, LanHEP prints it's 
prompt and waits for the keyboard input. In the last case, user's input is copied into the file lhep. log and 
can be inspected in the following. To finish the work with LanHEP, type 'quit. ' or simply press ~D (or ~Z 
at MS DOS computers). 

10.1 Options 

Possible options, which can be used in the command line to start LanHEP are: 
-OutDir directory Set the directory where output files will be placed. 

-InDir directory Set the default directory where to search files, which included by read and use state- 
ments. 

-tex LanHEP generates LaTeX files instead of CompHEP model tables. 

-f rc If -tex option is set, forces LanHEP to split 4-fermion and 4-color vertices just as it is made for 
CompHEP files. 

-texLines num Set number of lines in LaTeX tables to num. After the specified number of lines, LanHEP 
continues writing current table on the next page of LaTeX the output. Default value is 40. 

-texLineLength num Controls width of the Lagrangian table. Default value is 35, user can vary table 
width by changing this parameter. 

11 Default objects 

When LanHEP starts, it has already declared some frequently used symbols. They are: 

• specials gamma and gamma5 which are Dirac's 7-matrices; 

• special moment is replaced by the momentum of the corresponding particle; 

• let-substitution deriv is defined as -i*moment; 

• group color with indices types c3/c3b and c8; 
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specials lambda and f _SU3 which are Gell-Mann matrices and structure constant of SU(3) group; 

taul , tau2, tau3 are T-matrices t±, t%, T3; 

tau is a vector f = (ti, T2, T3); 

taup and taum are matrices = (r 1 ± ir 2 )/V2; 

taupm is a vector (t + ,t 3 ,t~); 

eps is the antisymmetrical tensor e (e 123 = 1); 

Taul, Tau2, Tau3 are generators of SU(2) group adjoint representation (3-dimensional analog of 
matrices) T 1 ,! 12 ,! 13 with commutative relations [Ti,Tj] = —ieiji-T}.: 

1 I -1 \ 1 I i \ / 10 °\ 

T 1 = — I — 1 1, T 2 = — [ -i -i , T 3 = ; 
^ 2 \0 10/ ^ 2 V0z0/ V 00 - 1 / 

Taup and Taum corresponds to T± = (T 1 ± iT 2 )/V2; 
Taupm is a vector T = (T + ,T 3 , T~). 
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12 Standard Model 

In this chapter we give an example of the LanHEP input file introducing the Lagrangian of the Standard Model 
in the t'Hooft-Feynman gauge. 

We start the description of the SM Lagrangian from the declaration of all parameters, see Fig. ||. The 
comments which follow each parameter explain its meaning. 



model 'Standard Model '/4. 
parameter 

EE = 0.31333: 'Elementary electric charge', 

GG = 1 . 117: 'Strong coupling constant (Z point) (PDG-94) ' , 

SW = 0.4740: 'sin of the Weinberg angle (PDG-94, "on-shell") ' , 

sl2 = 0.221: 'Parameter of C-K-M matrix (PDG-94)', 

s23 = 0.040: 'Parameter of C-K-M matrix (PDG-94)', 

sl3 = 0.0035: 'Parameter of C-K-M matrix (PDG-94)'. 

parameter CW = sqrt (1-SW**2) : 'cos of the Weinberg angle'. 

parameter 

cl2 = sqrt (l-sl2**2) : 'parameter of C-K-M matrix', 
c23 = sqrt (l-s23**2) : 'parameter of C-K-M matrix', 
cl3 = sqrt (l-sl3**2) : 'parameter of C-K-M matrix', 
parameter 

Vud = cl2*cl3: 'C-K-M matrix element', 
Vus = sl2*cl3: 'C-K-M matrix element', 
Vub = sl3: 'C-K-M matrix element', 

Vcd = (-Sl2*c23-cl2*s23*sl3) : 'C-K-M matrix element', 
Vcs = (cl2*c23-sl2*s23*sl3) : 'C-K-M matrix element', 
Vcb = s23*cl3: 'C-K-M matrix element', 
Vtd = (sl2*s23-cl2*c23*sl3) : 'C-K-M matrix element', 
Vts = (-cl2*s23-sl2*c23*sl3) : 'C-K-M matrix element', 
Vtb = c23*cl3: 'C-K-M matrix element'. 
0rthMatrix( {{Vud, Vus , Vub} , {Vcd, Vcs , Vcb} , {Vtd,Vts,Vtb}} ). 



Figure 3: Standard Model: parameters. 



Than we go to input particles of the model. These are vector gauge bosons, fermions which are leptons and 
quarks of three generations, and scalar Higgs boson (see Fig. |^). 

The next section of the input file (Fig. ||) introduces some useful let-substitutions. First, we enter doublets 
of leptons and quarks; for the quarks we supply also doublets rotated by CKM matrix. 

In the textbooks the physical gauge bosons W ± , Z, A are defined through the SU(2) gauge triplet W a (a = 
1,2,3) and U(l) field B as 

W ± = (W^iW 2 )/^, 
Z = W 3 cosO w - Bsiii0 w , 
A = W 3 sm9 w + Bcos9 w , 

where Ow is the Weinberg angle. If we would like to use the fields W a and B we have to declare them (as 
let-substitutions) through the physical fields W ± , Z and A by inverting these relations^. At last, we define EW 
charges g — e/ sin 6w and g\—ej cos#w. 

Now we can enter the Lagrangian terms. First, we define the EW gauge bosons Lagrangian (see Fig. ^) in 
the form 

where 

F^ = d M B v - d v B„ and = d^W? - d v W« - ge abc W*W, 
and for the gluon Lagrangian we use the form 

r QC D _ 1 fiaiiv r^a 
G - 5 



Note, that Bl stands for the B gauge field since B is already defined as b-antiquark. 
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vector 

A/A: (photon, gauge), 

Z/Z: ('Z boson', mass MZ = 91.187, width wZ = 2.502, gauge), 
G/G: (gluon, color c8, gauge), 

'W+'/'W-': C'W boson', mass MW = MZ*CW, width wW = 2.094, gauge). 

spinor 

nl : (neutrino , left) , el : (electron, mass Me = 0.000511), 

n2 : ( 'mu-neutrino ' , lef t) , e2:(muon, mass Mm = 0.1057), 

n3: ( 'tau-neutrino' ,left) , e3 : ( 'tau-lepton' , mass Mt = 1.777). 

spinor 



u 


( 


u- quark' 


, color 


c3), 








d 


( 


d- quark' 


, color 


c3), 








c 


( 


c- quark' 


, color 


c3, mass 


Mc = 


1 


300) , 


s 


( 


s- quark' 


, color 


c3, mass 


Ms = 





200) , 


t 


( 


t- quark' 


, color 


c3, mass 


Mtop 




170, width wtop = 1.442), 


b 


( 


b-quark' 


, color 


c3, mass 


Mb = 


4 


300 ). 



scalar H/H:(Higgs, mass MH = 200, width wH = 1.461). 



Figure 4: Standard Model: particles. 



let ll={nl,el}, L1={N1,E1}. 
let 12={n2,e2}, L2={N2,E2}. 
let 13={n3,e3}, L3={N3,E3}. 

let ql={u,d}, Q1={U,D}, qla={u, Vud*d+Vus*s+Vub*b} , Qla={U, Vud*D+Vus*S+Vub*B} . 
let q2={c,s}, Q2={C,S}, q2a={c,Vcd*d+Vcs*s+Vcb*b}, Q2a={C , Vcd*D+Vcs*S+Vcb*B} . 
let q3={t,b}, Q3={T,B}, q3a={t , Vtd*d+Vts*s+Vtb*b} , Q3a={T,Vtd*D+Vts*S+Vtb*B} . 
let Bl= -SW*Z+CW*A, W3=CW*Z+SW*A, Wl=( 'W+'+'W- ' ) /Sqrt2 , W2 = i* ( ' W+ ' - ' W- ' ) /Sqrt2 . 
let WW={'W+', W3, 'W-'}, WWc={'W-', W3, 'W+'}, WW1={W1, W2, W3} . 

let ghl=('W+.c'+'W-.c')/Sqrt2, gh2= i* ( ' W+ . c ' - ' W- . c ' )/Sqrt2 , gh3= CW*'Z.c'+SW*'A.c' , 
gh={ghl,gh2,gh3}. 

let Ghl = ('W+.C'+'W-.C')/Sqrt2, Gh2=i* ( ' W+ . C ' - ' W- . C ' ) /Sqrt2 , Gh3= CW*'Z.C'+SW*'A.C , 

Gh={Ghl,Gh2,Gh3}. 
let g=EE/SW, gl=EE/CW. 



Figure 5: Standard Model: useful substitutions. 

where 

is the gluon field, g s is a strong interaction coupling. 



let tB~mu~nu=deriv~mu*Bl~nu-deriv~nu*Bl~mu . 

let tW~mu~nu~a=deriv~mu*WWl~nu~a - deriv~nu*WWl~mu~a - g*eps"a"b"c*WWl~mu~b*WWl"nu~c . 
let tG~mu~nu~a=deriv~mu*G~nu~a - deriv~nu*G~mu~a - i*GG*f ~a~b~c*G~mu~b*G~nu~c . 
lterm -tB**2/4-tW**2/4-tG**2/4 . 



Figure 6: Standard Model: gauge field Lagrangians 

The next step is the interaction of fermions with gauge bosons (Fig. [?]). First, we define interaction of 
left-handed fermions with W a and B fields, taking into account CKM mixing of quarks; the corresponding term 
is 

h^-^f- fa - \r a W; - gtYB^ ^ 

where tp are left-handed lepton and rotated by CKM matrix quark doublets, Y is a hypercharge. Than we 
declare the right-hand fermions interaction with B field; in this case we can take quarks without CKM mixing, 
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since the interaction of quarks with Z-boson and photon is flavor-diagonal. Lagrangian term is 



where if) are right-handed lepton and quark singlets. 



lterm anti (psi) *gamma* ( l-g5) /2* (i*deriv- 


-g*taupm*WW/2-Y*gl*Bl) *psi 


where 


psi=ll, Y= -1/2; 






psi=12, Y= -1/2; 






psi=13, Y= -1/2; 






psi=ql, Y=l/6; 






psi=q2, Y=l/6; 






psi=q3, Y=l/6. 






lterm anti (psi) *gamma* (l+g5) /2* (i*deriv- 


-Y*gl*Bl)*psi 




where 






psi=el, Y= -1; 






psi=e2, Y= -1; 






psi=e3, Y= -1; 






psi=u, Y=2/3; 






psi=c, Y=2/3; 






psi=t, Y=2/3; 






psi=d, Y=-l/3; 






psi=s, Y=-l/3; 






psi=b, Y=-l/3. 







Figure 7: Standard Model: gauge interaction of fermions. 

The next step is the introduction of the Higgs sector: interaction of Higgs doublet expressed through physical 
field H and unphysical goldstone fields Zf, Wj as 



$ = 



-iWf 



-L(v + H + iZ f ) 



where vacuum expectation value is v = 2M\y/g- Unphysical goldstone fields read in LanHEP as 'W+.f ' , 
'W-.f ' , 'Z.f ' . In Fig. H we enter self-interaction terms and interaction of this field with gauge bosons as: 

Lm aa s = (D^)*(D^) + 2A($*$ - v 2 /2) 2 , 
where D„ = + zf T a W£ + ifB^ X = 



let phi = { -i*'W+.f, (vev(2*MW/EE*SW)+H+i* , Z.f ')/Sqrt2 }, 
Phi = { i*'W-.f, (vev(2*MW/EE*SW)+H-i*'Z.f ')/Sqrt2 }. 

lterm -2*lambda* (pp*PP+v**2/2) **2 

where lambda=(EE*MH/(MW*SW))**2/16, v=2*MW*SW/EE . 

lterm (DPhi) * (Dphi) where 

DPhi = (deriv~mu-i*gl/2*Bl~mu) *Phi~a - i*g/2*taupm~a~b~c*WWc~mu~c*Phi~b, 
Dphi = = (deriv~mu+i*gl/2*Bl~mu) *phi~a + i*g/2*taupm~a~b~c*WWc~mu~c*phi~b . 



Figure 8: Standard Model: self-interaction of Higgs field and its interaction with gauge bosons. 
Yukawa terms are shown in Fig. |^. First, we generate lepton masses with terms 



V2AL 



- 1 + T 5 - 1 - 7 5 

ip a ' -i/"£ a + i> -if> a <!>* a 



where ip a is lepton doublet, if) is down lepton and Mi is its mass. The term is applied to the leptons of second 
and third generations (electron is massless in this model). 
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Next we do the same thing with down quarks. Here we must take into account the CKM mixing, so the 
previous formula takes the form 

n gM l J /_l + 7 5 -I-7 5 

Ly d = ~ " I W^LJ-^ + ^.__L^«$*« 



where M*- 7 = diag(Md, M s , M^) U 3k , U being the CKM matrix; ipf and are quark doublets and down 
singlets, indices i,j numerate generations. 

At last we generate masses of upper quarks with the terms 

gM q {Ta 1 + 7 5 .ab, ifl *b , J. 1 -^ J.a.ab^b 



for three generations of quarks with M q = (M u , M c , M t ) 



lterm -M/MW/Sqrt2*g* (anti (pi) * (l+g5) /2*pr*phi + anti(pr)*(l-g5)/2*pl*Phi ) where 
M=Mm, pl=12, pr=e2; 
M=Mt , pl=13, pr=e3. 

lterm -M/MW/Sqrt2*g* (anti (pi) * (l+g5) /2*pr*phi + anti(pr)*(l-g5)/2*pl*Phi ) where 

M=Vud*0, pl=qla, pr=d; °L here and below stands for Md 

M=Vus*Ms, pl=qla, pr=s; 

M=Vub*Mb, pl=qla, pr=b; 

M=Vcd*0, pl=q2a, pr=d; 

M=Vcs*Ms, pl=q2a, pr=s; 

M=Vcb*Mb, pl=q2a, pr=b; 

M=Vtd*0, pl=q3a, pr=d; 

M=Vts*Ms, pl=q3a, pr=s; 

M=Vtb*Mb, pl=q3a, pr=b. 
lterm -M/MW/Sqrt2*g*(anti(pl)*(l+g5)/2*i*tau2*pr*Phi + anti (pr) * (l-g5) /2*i*pl*tau2*phi ) 
where 

M=0, pl=qla, pr=u; 

M=Mc, pl=q2a, pr=c; 

M=Mtop, pl=q3a, pr=t . 



Figure 9: Standard Model: Yukawa terms. 

The final stage is the introduction of gauge fixing terms and terms with Faddev-Popov ghosts. 
Gauge fixing terms in the t'Hooft-Feynman gauge read as 

Lgf = -\{d^) 2 - i(9 M G£) 2 - + MwWfWuW-" + M W WJ) - \{d^ + M z Z f f, 



the corresponding LanHEP code is shown in Fig. 10 



lterm -1/2* (deriv*A) **2 . 
lterm -1/2* (deriv*G) **2 . 

lterm -1/2* (2* (deriv* »W+' +MW* 'W+ . f ' ) * (deriv* 'W- ' +MW* 'W- . f ' ) + 
(deriv*Z+MW/CW* ' Z . f ' ) **2) . 



Figure 10: Standard Model: gauge fixing terms. 
The interaction of ghost fields reads as (we omit here bilinear terms, since they are not treated by CompHEP): 

L G FP = ig s f ahc C G G b ^C G 

for gluons, and 

l fp =-9£ *^ w W fi C w 
for the interaction of EW ghosts with gauge bosons, and finally 
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+HC Z C Z / cos 2 6 W ~ 2iZ f C+C-) + 

+ o „ (Wf((l - 2 sin 2 e w )C-C z + C^Cz + 2 sin6 w cos W C-C A ) ~ 

2 cos Ow sin 6% J 

-W7((l - 2 sin 2 6W)C+C*z + C + (7 Z + 2sin0 w a»0 w -<7+C , . A )), 

Here (a=1..8) are ghost fields corresponding to gluon, C w (a=1..3) are ghost fields of SU{2) gauge 
bosons, Cb corresponds to U(l) gauge field and 

Cz = C w cos 9w — Cb sin 9 W , 
Ca — C w sin Gw + Cb cos 9 W . 

The corresponding part of LanHEP input file is shown in Fig. HJ. 



lterm i*GG*f *ccghost (G) *G*deriv*ghost (G) . 
lterm i*g*eps*Gh*WW*deriv*gli. 

lterm -1/2* (2* (deriv* >W+' +MW* 'W+ . f ' ) * (deriv* 'W- ' +MW* 'W- . f ' ) + 

(deriv*Z+MW/CW* ' Z .f ' ) **2) . 
lterm -MW*EE/2/SW* ( (H+i* ' Z . f ')*('W-.C'*'W+.c' + ' W+ . C ' * ' W- . c ' ) 

+H* ' Z . C ' * ' Z . c ' /CW**2-2*i* ' Z . f ' * > W+. C ' * ' W- . c ' ) . lterm 
i*EE*MW/2/CW/SW*('W+.f ' * ( 'W- .C' * ' Z . c ' * (1-2*SW**2)+'W- . c ' * ' Z . C' 
+2*CW*SW*'W-.C'*'A.c') - 
'W-.f '*('W+.C'*'Z.c , ^l-2*SW**2) + , W+.c'*'Z.C' 
+2*CW*SW* > W+ . C ' * ' A . c ' ) ) • 



Figure 11: Standard Model: terms with Faddev-Popov ghost fields. 

To complete the description of Standard Model, we prescribe also LaTeX names for particles and parameters, 
see Fig. |l2] . 



SetTexName ( [u=u , U= ' \\bar{u} ' , d=d, D= ' \\bar{d} ' ] ) . 

SetTexName ( [s=s, S='\\bar{s}', b=b, B='\\barb']). 

SetTexName ( [c=c, C='\\bar{c}', t=t, T=' \\bar{t}'] ) . 

SetTexName ( [el=e, El=' \\bar{e} ' , nl='\\nu~e', Nl='\\bar{\\mi}~e'] ) . 

SetTexName ( [e2='\\mu' , E2=' \\bar{\\mu}' , n2='\\nu"\\mu' , N2='\\bar{\\nu}~\\mu'] ) . 

SetTexName ( [e3='\\tau' , E3= ' \\bar{\\tau} ' , n3=' \\nu~\\tau' , N3= ' \\bar{\\nu}"\\tau' ] ) . 

SetTexName ( [EE=e , GG='g_s', SW='s_w', CW='c_w', MZ='M_Z']). 

SetTexName ( [Me='M_e' , Mm='M_\\mu' , Mt='M_\\tau>, Ms='M_s', Mc='M_c', Mb=M_b, Mtop=M_t] ) . 



Figure 12: Standard Model: LaTeX names for particles and parameters. 
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Appendix 



A Processing 4-color vertices 

CompHEP Lagrangian tables don't describe explicitly color structure of a vertex. If color particles present in 
the vertex, the following implicit convolutions are assumed (supposing p, q, r are color indices of particles in the 
vertex) : 

• S pq for two color particles A p , A 2 ; 

• Xp q for three particles, which are color triplet, antitriplet and octet; 

• f pqr for three color octets. 

Other color structures are forbidden in CompHEP. 

So, to introduce the 4-gluon vertex f pqt 'G < f 1 G^f pst G s l]j G t v one should split this 4-legs vertex into 3-legs vertices 

fpqrnqnr VP ■ 




Here the field X? is Lorentz tensor and color octet, and this field also has constant propagator. If gluon 
name in CompHEP is 'G', one can use name 'G.t' for this tensor particle; its indices denoted as 'm_' and 
'M_' ('_' is the number of the particle in table item). 

B Feynman rules for the Standard Model 

Here we present the tables generated by LanHEP in LaTeX format for the Standard Model (see Section 12). 



Table 2: Standard model: particles. 



p 


aP 


Name 


Spin 


EM charge 


Color 


2-lcgs vertex 


Comment 




A v 


photon 


1 





1 




gauge 




z v 


Z boson 


1 





1 


-?(6irfr-w 2 .r) 


gauge 


G m 


Gyq 


gluon 


1 





8 


-pfp^Spq 


gauge 




w- v 


W boson 


1 


1 


1 


-q pv (pfpf - MW 2 ) 


gauge 




v e b 


neutrino 


1/2 





1 


-M (lab 7£c7 c 5 6 ) 


left 


e a 


eb 


electron 


1/2 


-1 


1 


IX u 
-Vxlab 




" a 


y"b 


mu-neutrino 


1/2 





1 


-M Kb-^acllb) 
-(P^ab + M^Sab) 


left 


[i.a 


Ub 


muon 


1/2 


-1 


1 




V T a 


v T b 


tau-neutrino 


1/2 





1 


-iPl (lib-lLllb) 
-(P^b + Mr-Sab) 


left 


T a 


fb 


tau-lcpton 


1/2 


-1 


1 




U a p 


U_bq 


u-quark 


1/2 


2 

_! 


3 


-Pl 6 pqlab 




dap 


dbq 


d-quark 


1/2 


3 


-\pt5 pq 2^ b 




Cap 


Cbq 


c-quark 


1/2 


_! 


3 


Spq(P^ b + M c -6 ab ) 




s ap 


Sbq 


s- quark 


1/2 


3 


-±S pq (2p^ b + 2M s -5 ab ) 
(Pilab + M t ■ Sab) 




tap 


tpq 


t-quark 


1/2 


I 

3 


3 




bap 


hq 


b-quark 


1/2 


3 


-±S pq (2p^ b + 2M b -5 ab ) 
-(MH 2 -pW) 




H 


H 


Higgs 








1 
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Table 3: Standard Model: parameters. 



Parameter 



Value 



Comment 



EE 

GG 

SW 

sl2 

s23 

sl3 

CW 

cl2 

c23 

cl3 

Vud 

Vus 

Vub 

Vcd 

Vcs 

Vcb 

Vtd 

Vts 

Vtb 

MZ 

wZ 

MW 

wW 

Me 

Mm 

Mt 

Mc 

Ms 

Mtop 

wtop 

Mb 

MH 

wH 



0.31333 

1.117 

0.474 

0.221 

0.04 

0.0035 

Sqrt(l-SW**2) 

Sqrt(l-sl2**2) 

Sqrt(l-s23**2) 

Sqrt(l-sl3**2) 

cl2*cl3 

sl2*cl3 

sl3 

-sl2*c23-cl2*s23*sl3 

cl2*c23-sl2*s23*sl3 

s23*cl3 

sl2*s23-cl2*c23*sl3 

-cl2*s23-sl2*c23*sl3 

c23*cl3 

91.187 

2.502 

MZ*CW 

2.094 

0.000511 

0.1057 

1.777 

1.3 

0.2 

170 

1.442 

4.3 

200 

1.461 



Electromagnetic coupling constant (1/128) 

Strong coupling constant (Z point) (PDG-94) 

sin of the Weinberg angle (PDG-94," on-shell"i 

Parameter of C-K-M matrix (PDG-94) 

Parameter of C-K-M matrix (PDG-94) 

Parameter of C-K-M matrix (PDG-94) 

cos of the Weinberg angle 

parameter of C-K-M matrix 

parameter of C-K-M matrix 

parameter of C-K-M matrix 

C-K-M matrix element 

C-K-M matrix element 

C-K-M matrix element 

C-K-M matrix element 

C-K-M matrix element 

C-K-M matrix element 

C-K-M matrix element 

C-K-M matrix element 

C-K-M matrix element 

mass of Z boson 

width of Z boson 

mass of W boson 

width of W boson 

mass of electron 

mass of muon 

mass of tau-lepton 

mass of c-quark 

mass of s-quark 

mass of t-quark 

width of t-quark 

mass of b-quark 

mass of Higgs 

width of Higgs 
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Fields in 


the vertex 


Variational derivative of Lagrangian by fields 






G pr 


o. f„„ r -(voa pp - v P .a pv - v*ia vp + v p a tJ,v + vHa" p - vKa pp \ 

ysJpqr\F3y F2a ^3 a ~ yiiJ ~ t"2. i) PlU J 






G p r G(j S 


gH^g^fputfrst - g^g up f Pq tfrst + g pv g pa f P rtf qst 

>y bf Jpstjqrt y y Jprtjqst y y J pstj qrtj 


w+p 


W~jj 


^p 


Suj \Piy i j \y Piy ~ f2y ~ysy p^y j 


Ap 




p 


-e(p p q^ u - p p q vp - p»q pp + p p q vp + p»q pp - p p q pv ) 


w+p 


w- v 




-^l(2q^ Q P' T - q pp q va - q pa q vp \ 
s -z \^y y y y y y j 


Ap 




VV p ZJ(j 


- c " g2 (2o pa o up — a pv a pa - a pp a va \ 

Sw v " y y y y y j 


Ap 


A v 


W + p W~ a 


v " y y y y y } 


w+p 




VV p VV (j 


£(2n pv n pcr — a pa n up — n pp n v(T \ 

\^y y y y y y ) 






z„ 

p 


4 c w s w \ ' ab lac lebj 


u a 




w+ , 

w p 


4 s w \ <ab lac 1 cb I 


e a 


7 ,c 
V b 


w~ 

VV p 


4 s w \ 'ab lac 1 cb I 


e a 


e-b 


z„ 

p 


1 6 ( 4:S 2 ^ P +^ P ^ P l\) 

4 c m s m V w lab ~ lab lac IcbJ 


e a 


eb 


A 

P 


M 

t 'lab 


CM 
u a 


V^b 


z„ 


4 c^s^ V 'ao lac IcbJ 


HP 

u a 


Mb 


w+, 

vv p 


4 s m V lab lac 1 cb 1 


Ma 


y^b 


vv p 


— L / P _ p 5 ) 

4 s m V 'ab /ac /cb/ 


Ma 


Mb 


z 

p 


4 c„s„ V lab ~ lab lac leb) 


Ma 


Mb 


A 

P 


M 

t 'lab 


V T a 


v T b 


z„ 


1 e / M 7 (i 7 M 
4 c w s w \ 1 ab lac IcbJ 


V T a 


Tb 


w+„ 

w p 


4 s m V 'ab 'ac /cb/ 


T a 


v T b 


vv p 


4 s m V 'ab lac IcbJ 


f a 


Tb 


z 

p 


4c„s„V ^s w J ab -tJ ab laclcb) 


fa 


Tb 


An 

P 


M 

c 'ab 


Hap 


Ubq 


z, 

p 


12 c„s„ °P9 V °*w 'ab ^ J 'ab ° 'ac 'cb/ 


Hap 


Ubq 


4 

M 


--e<5 t m 

3 cu pq lab 


^ap 


dbq 


W+ 
" M 


1 ev^V«d r / M u 5 \ 
4 s„ °P9 V 'ab 'ac /cb/ 


Hap 


Sbq 


W+ 

M 


4 s„ °P° V 'ab lac IcbJ 


Hap 


hq 




leV2Vub 5 ( H V»~5\ 
4 s„ "P9 V 'ab 'ac 'cb/ 


dap 


Ubq 


" P 


4 s„ V 'ab 'ac 'cb/ 


d a p 


dbq 


z, 

p 


12 c„s ra P9V w lab ' ° lab ° lac IcbJ 


dap 


dbq 


A 

P 


-eS 2-y p 

6 e °P° 'ab 


Sap 


Ubq 




leV2Vus S J P ^P,y5) 
4 s w "Ply lab lac IcbJ 


$ap 


Sbq 


z„ 

p 


1 e 6„„ ( 4s 2 -7^ + 3-7^ 3t^ •y\) 

12 c„s ra P°V w 'ab ' ° lab ° 'ac 'cb/ 


Sap 


Sbq 


,4 


-eS 27 M 

6 e0 P° Z 'ab 


&ap 


Ubq 


" M 


lev^y«br / p „,p 5) 
4 s w PQ V 'ab lac IcbJ 


&ap 


bbq 


L P 


12 c,„s,„ 4S «>7ah + ^7ah J 7ac7 c bJ 



t lelds m 


the vertex 


Variational derivative of Lagrangian terms by fields 


bap 


hq 


A„ 
n 






Cap 


Cbq 


Z„ 
n 


-T2^p q (-^ab + 'Kb 


- ZiLilb) 


Cap 


Cbq 


A, 


-^efipqlab 




Cap 


dbq 


w + , 

yv P 


1 eV2Vcd g In a 5 \ 
4 s„ °Pq\lab laclcb) 




Cap 


Sbq 


yy n 


leV2Vcs S I n u 5) 
4 s„ U PQ \ lab lac leb) 




Cap 


bbq 


yy p 


leV2Vcbx IP u 5\ 
4 s w °Pq\lab laclcb) 




dap 


Cbq 


y y 


1 eV2Vcdz I n u „5 \ 
4 s w °P1\lab laclcb) 




$ap 


Cbq 


yy fj. 


4 s w °P1\lab laclcb) 




&ap 


Cbq 


yy fj. 


-Y-^^qUb-lLllb) 




tap 


tbq 


Z, 

n 




- ZlZcllb) 


tap 


tbq 


A, 


— ^eSpqla b 




tap 


dbq 


W + „ 
yy p 


i eV2Vtd s i n n ~5 \ 

4 s m U P° \ lab lac ) cb) 




tap 


Sbq 


W + , 
yy p 


icV2Vts 5 in y^V) 

4 s„ U P1 \ >ab lac 1 cb) 




tap 


hq 


W+, 
yy p 


1 cV2Vtb s i n .,5 \ 
4 s w °P1 V lab laclcb) 




dap 


tbq 


yy n 


1 cV2Vtd § i n n 5) 
4 s m °P<l\lab laclcb) 




Sap 


tbq 


yy n 


leV2Vts 6 i n ~n^5\ 
4 s„ U P1 V lab lac leb) 




bap 


tbq 


yy n 


lcV2Vtb S i n ^.5\ 
4 s w u P1\lab laclcb) 




U a p 


Ubq 


G [x r 


9slab X pq 




d a p 


dbq 




9slabKq 




Cap 


Cbq 


G p,r 


9slab X pq 




Sap 


Sbq 


G p,T 


9slab X pq 




tap 


tbq 


G pr 


9slab X pq 




bap 


hq 


G p,r 


9slab X pq 




U a p 


Sbq 


w + 

yy f 


4 Mffs, U Pq \ Uab i lab) 




Sap 


Sbq 


H 


1 eM 3 x X 

2 Mffs, °abOpq 




Sap 


Sbq 


Zt? 

ZJ 


1 ieM s x 5 

2 Mffs„ u Pq lab 




Sap 


Ubq 


w~ 

yy f 


1 ieM s V2Vus x (X „5 \ 




Uap 


hq 


wi 

yy f 


1 zcM b V2Vubx (X , „,5 \ 
4 Mffs, Vpq\Vab + 7 ah J 




bap 


Kq 


H 


2 MWs„ u ab u pq 




bap 


hq 


Zv 


1 ieM b r ,,5 

2 MWs„ lab 




bap 


Ubq 


yy p 


1 jeM,V^V«6 e /c 5 \ 
4 MlVs ra Vpq\Vab lab) 




Cap 


Sbq 


F 




M c 8 ab + M cl 5 ab ) 


Sap 


Cbq 


yy p 


-IlMf S pq( M sSab-M s ^ ab 


- M c S ab - M cl 5 ab ) 


Cap 


hq 


yy p 


\lM^ S pq( M ^ + M bl 5 ab - 


M c S ab + M cl l b ) 


bap 


Cbq 


yy p 




- M c S ab - M c7a \) 


tap 


Sbq 


W+ 


\ li T^Spq{Ms8ab + M sl l b - 


M t S ab + M tl l b ) 
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fields m 


the vertex 


Variational derivative of Lagrangian by fields 


_ 

$ap 


v oq 


w F 




lie ,y^y tS S nn (MJnh M.7 5 , MfSnh Mt^u) 

4 MWs m PI \ s a0 s lab lvJ -t u ab lrJ -t l a i,) 


T 

v ap 


hi,, 
u bq 


w+ 




1 ieV2Vtb r J M r , M 5 Mt8„h + M+tM 


h 

u ap 


^bq 


Wp 




4 Mffs, °P9V w &°ah m bl ab lvl t O a b lvl tl a b) 




Cbq 


H 




1 eM c r c 

2 MWs m UabU pq 


^ap 


^bq 


Z F 




1 ieM c r 5 

2 MWs w V P1 'ab 


^ap 


^bq 


Wp 




lieM c V2Vcd 5 ( 5 , 5\ 
4 MWs m u pq\ u ab^ l a b) 


7> 

^ap 


Ubq 


w+ 




1 ieM c \/2Vc<i r /r 5 \ 
4 MWs m u pq\ u ab J a b) 


^ap 


^bq 


if 




1 eMt S nh 5 nn 


f 

^ap 


ft, 
°bq 


Z F 




1 ieM t r 5 

2 Mffs, °P a 'ab 


d 

Map 


tu 
L bq 


Wp 




lieM t V2Vtd g ( S , +ry 5\ 
4 MWs w °pq\°ab "t" / a hj 


^ap 


(thn 


w+ 




4 MWsu, u pq\ u ab lab) 


v a 




w+ 




lieM„V2( 5 . 5 \ 
4 MVKs„ V 0ab T 'afcj 


77 

ra 


H'b 


if 




1 elf, ^ 


Uj„ 
h^a 


Uh 
f-*0 


Z F 




1 "M, 5 

2 Mffs. 'aft 


ra 


1 


Wp 




4 Mffs, V Uab 'ab/ 


V a 


Tb 


w+ 




1 ieM T ^/i , 5 \ 
4 Mffs, V" ab ' lab) 


1 a 


' b 


if 




1 eM T e 

2 MWs w Uab 


' a 


Th 


Z F 




1 icM T 5 

2 Mffs. 'ab 


' a 


1 


Wp 




4 Mffs. V Uab 'ab/ 




wt 

rv p 


Wp 


w F 


1 e 2 MH 2 

2 MW^sJ, 


if 


wt 

VV p 


Wp 




1 eMH ! 

2 Mffs. 


if 


if 


w+ 


w F 


1 e 2 MH 2 
4 MW 2 4 


Wt 


w„ 

VV p 


Z F 


Z F 


1 e 2 MH 2 

4 Miy 2 4 


if 


H 


if 




3 eMH 2 
2 MWs w 




/j t 


Z F 




1 eMH 2 

2 Mffs. 


if 


H 


if 


if 


3 e 2 MH 2 

4 MW 2 4 


if 


if 


Z F 


Z F 


1 e 2 MH 2 
4 MW 2 sJ, 




b 


Z f 


Z f 


3 e 2 MH 2 

4 MH^sJ, 


wt 

VV p 


wz 

VV p 


z. 




1 6 (2 s 2 ?/ 1 B M 2S 2 B M +B A ' N | 


if 


z 


Z F 




2 c w s w \ -TPl) 


A 

M 


W + 
VV p 


w F 




c \P3 P2 ) 


if 


w + 

VV p 


w-» 




2 s ra V^l ^2 ) 




w~ 

VV y, 


Z F 




2 s ra V^3 ^1 j 


if 


w + , 

VV fl 


w F 




2 s„ V^3 ^1 / 




ll p 


Z F 




2 is™ V^2 -P3 / 


w + 


w~ 

VV p 




z v 


1 f 2 , o^^^s 4 4s 2 4- 
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£ lelds m 


the vertex 


Variational derivative of Lagrangian by fields 


if 




z v 




eMW uv 




H 




z v 


1 2 

2 c^s^ y 




z v 


Z F 


Z F 


2 c^s^ y 




w+ 


w F 


z v 


e2 a'"' (2 s 2 \ 

cn,s w y \ w > 




w~, 


z v 




ieMW s w iiv 


If 


w+ 


w-. 


z v 


1 i£_ n V 

2 c w V 


VK+ 


w~. 


z v 


Z F 


2 Cw y 




w F 


z v 




ieMWs w iiu 


if 




w F 


z v 


2 c m f 




w F 


z v 


Z F 


2 c m y 




A, 


w+ 


w F 


2p 2 o^ u 




w+ 


w- v 




-ieMWgV v 


A, 


H 


w+ 


w- v 


-h — Q^ 

2 s%y 


A, 


w+ 


w- v 


Z F 


1 2 

— ke_a UV 

2 s w y 


A, 


w+ u 


w F 




ieMWgV 


A, 


H 




w F 


i ■ 2 

1 le iiv 

2 s„ y 


A, 




w F 


Z F 




w+ u 


w+ 


w- v 


w F 


i 2 

2^T 5 


H 




w- v 




eMW uv 

s w y 


H 


H 


w+» 


w- v 


2 


w+ u 


w- v 


Z F 


Z F 


i 2 

1 e fjtv 




c G q 


c^~* r 






c w+ 


w-. 


c z 




« P3 


C A 


c w+ 


w-^ 




e Pl 


c w+ 


C w- 


z. 




c m e u 


A, 


c w+ 


c w - 




e P3 


w+ u 


c w- 


c z 




_c M e V 


C A 


w+ u 


c w- 




pn^ 

—ep 1 


c w + 


c w- 


z» 






Aft 


(jW+ 


c w ~ 




e P2 


w+ u 


C w- 


c z 






c w + 




c z 






c z 


w+ u 


c w ~ 




prf 
e P3 


c z 




w-^ 




M 

~ e P2 


H 


(jW+ 


c w - 




1 eMW 

2 s„ 


H 


C w + 


c w ~ 




1 eMW 

2 s 1( , 
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Fields in the vertex 


Variational derivative of Lagrangian by fields 




c w- 


Z F 


1 ieM W 

2 s„ 




c w- 


Z f 


1 ieMW 

2 s w 


H 


C z 


C z 


1 eMW 

2 cj,s„ 


w+ 




c z 


1 ieMW ( o„2\ 

2 c w s w \ Ab w) 


w+ 


C w ~ 


c z 


1 ieMW 

2 c„s ro 


c z 


w+ 


C w- 


leivi vv 


(jW+ 


w- 


C z 


1 ieMW I o„2\ 

2 c„s„ V Z,b w) 


(jW+ 


Wp 


c z 


1 ieMW 

2 c m s m 


C z 


(jW+ 


Wp 


—ieM W 
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